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ABSTRACT 

We report the discovery from WASP-South of a new Jupiter-like extrasolar planet, WASP- 
16b, which transits its solar analog host star every 3.12 days. Analysis of the transit photometry 
and radial velocity spectroscopic data leads to a planet with i?, p = 1.008 ± 0.071 Rj up and M p = 
0.855±0.059 Mj up , orbiting a host star with R* = 0.946 ± 0.054 R Q and M* = 1.022±0.101 M . 
Comparison of the high resolution stellar spectrum with synthetic spectra and stellar evolution 
models indicates the host star is a near-solar metallicity ([Fe/H]= 0.01 ± 0.10) solar analog 
(T cff = 5700 ± 150 K, log g= 4.5 ± 0.2) of intermediate age (r = 2.3±|;| Gyr). 



Subject headings: planetary systems : individual: WASP-16b — stars: individual () — stars: abundances 
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1. Introduction 

There are currently over 300 known exoplan- 
eta^l with the majority of them discovered through 
the radial velocity technique. A growing number 
of exoplanets in recent years have been discovered 
through the transit method. Transiting exoplan- 
ets are particularly valuable as they allow param- 
eters such as the mass, radius and density to be 
accurately determined and further studies such as 
transmission spectroscopy, secondary eclipse mea- 
surements and transit timing variations to be car- 
ried out. 

There are several wide angle surveys that have 
been successful in finding transiting exoplanets 
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aroun d brig ht stars, namely H AT ( Bakos et al 



20021) . TrES (lAlonso et alJlioO 



McCullough et al 



2005) and WASP (IPollacco et all |2006j). The 



WASP Consortium conducts the only exoplanet 
search currently operating in both hemispheres 
although HATnet is planning a southern exten- 
sion and several groups are planning searches 



Strassmeier et al. 2007 



from Antarctica (e . g . 
Crouzet et aDl2009h . 

We report the discovery from the WASP-South 
observatory of a ~ 0.86 Mj up mass companion or- 
biting a V ~ 11.3 close solar analog WASP-16 
(=TYC 6147-229-1, USNO-B1.0 0697-0298329). 

2. Observations 

2.1. Photometric observations 

WASP-South, located at SAAO, South Africa, 
is one of two Super WASP instruments and com- 
prises eight cameras on a robotic mount. Each 
camera consists of a Canon 200mm f/1.8 lens 
with an Andor 2048 x 2048 e2v CCD camera giv- 
ing a field of view of 7.8° x 7.8° and a pixel 
scale of 13.7". Exposure times were 30 s and the 
same field is returned to and reimaged every 8- 
10 minutes. Further details of the instrument, 
survey and da t a red uction pipelines are given in 
Pollacco et al. ( 20061 ) and the candidate selection 



proced ure i s described in ICollier Cameron et al . 
(|2007l ) and IPollacco et all (|2008l ) and references 
therein. 

WASP-16 was observed for a partial season in 
2006, a full season in 2007 and a further par- 
tial season in 2008 with the distribution of data 
points as 3324 points (2006), 6013 (2007) and 4084 
(2008). The 2007 light curve revealed the pres- 
ence of a ~ 1.3% dip with a period of ~ 3.11 days. 
The transit coverage in the other two seasons was 
very sparse, particularly in 2006, and there is only 
evidence for 2 partial transits in the 2008 data. 
WASP-16 was a fairly strong candidate for follow- 
up despite the sma ll number of transits, passing 
the filtering tests of lCollier Cameron et alJ ( 2006h 
with a signal to red noise ratio, SVcd = 9.38 (with 
SVcd > 5 required for selection), 'transit to anti- 
transit ratio' Ax 2 / Ax 2 - = 2.5 (Ax 2 / Ax 2 - > 1.5 
required for selection) and no measurable ellip- 
soidal variation. 

The SuperWASP light curve showing a zoom 
of the transit region, along with the model tran- 



sit fit, is shown in Figure [T] In order to better 
constrain the transit parameters, follow-up high 
precision photometric observations with the Swiss 
1.2m+EULERCAM on La Silla, were obtained in 
the I c band on the night of 2008 May 04 and are 
shown in Figure [2j 

2.2. Spectroscopic observations 

In order to confirm the planetary nature of the 
transit signal, we obtained follow-up spectroscopic 
observation with the Swiss 1.2m+CORALIE spec- 
trograph. The data were processed through the 
standard CORALIE redu c tion pipeline as de- 
scribed by iBaranne et al. ( 19961 ) with an addi- 
tional correction for the blaze function. Fourteen 
radial velocity measurements were made between 
2008 March 10 and 2008 August 04 and an ad- 
ditional sixteen between 2009 Feb 19 and 2009 
June 03 (see Table [T]) by cross-correlating with 
a G2 template mask. The resulting radial veloc- 
ity (RV) curve is shown in Figure [3] The low 
amplitude RV variation clearly supports the exis- 
tence of a planetary mass companion. In order to 
rule out a non-planetary explanation for the ra- 
dial velocity variation such as a blended eclipsing 
binary or starspots, we examined the line-bisector 
spans. Contamination from an unresolved eclips- 
ing binary will cause asymmetries in the spec- 
tral line profiles and line bisector span variations 



(jQueloz et alj|200lt [Torres et alj|2005l) . As can 



be seen from the lower panel of Figure [3l there 
is no sign of variation with phase of the bisector 
spans and their amplitude is much smaller than 
the radial velocity variation. This supports the 
conclusion that the radial velocity variations are 
due to a planet orbiting the star and not some 
other cause. 

3. WASP-16 System Parameters 
3.1. Stellar Parameters 

The individual CORALIE spectra are of rela- 
tively low signal-to-noise, but when co-added into 
O.OlA steps they give a S/N of around 70:1 which 
is suitable for a photospheric analysis of the host 
star. In addition, a single HARPS spectrum was 
used to complement the CORALIE analysis, but 
this spectrum had relatively modest S /N of around 
50:1. The standard CORALIE/HARPS pipeline 
reduction products were used in the analysis. 
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Table 1 

coralie radial velocities for wasp-16. 



Time of obs. Rad. Vel. orv Bisector span 
(BJD-2450000) (kms^ 1 ) (kms" 1 ) (kms" 1 ) 



4535. 


.864842 


-1 


.99772 


0. 


.01591 


0.00306 


4537, 


.849158 


-1 


.96688 


0. 


,00853 


-0.04553 


4538. 


.858364 


-2 


.00734 


0. 


,00899 


-0.03129 


4558. 


.780835 


-1 


.83336 


0. 


,00723 


-0.02779 


4560. 


.709473 


-2 


.00513 


0. 


,00725 


-0.02403 


4561. 


,688137 


-1 


.82730 


0. 


,00785 


-0.03998 


4589. 


.705102 


-1 


.84255 


0. 


,00875 


-0.04520 


4591. 


.706755 


-2 


.03571 


0. 


,00892 


-0.03221 


4652, 


.495906 


-1 


.82493 


0. 


,00808 


-0.03209 


4656, 


.551645 


-2 


.02421 


0. 


,00787 


-0.02555 


4657, 


.577293 


-1 


.96640 


0. 


.00957 


-0.01827 


4663 


.539741 


-2 


.02961 


0. 


.00969 


-0.02661 


4664, 


.616769 


-1 


.78590 


0, 


.01108 


-0.04350 


4682, 


.521501 


-1 


.98118 


0, 


.00754 


-0.02123 


4881, 


.869213 


-2 


.02245 


0, 


.00813 


-0.02760 


4882. 


.801025 


-1 


.83289 


0, 


.00823 


-0.03739 


4884. 


.737094 


-2 


.04565 


0, 


.00778 


-0.01672 


4891. 


.805707 


-1 


.90043 


0. 


.00798 


-0.01009 


4892. 


.723980 


-1 


.83413 


0. 


.00891 


-0.02116 


4941 


.728231 


-1 


.88737 


0, 


.00748 


-0.04134 


4943 


.730102 


-2 


.04677 


0, 


.00753 


-0.01825 


4944. 


.739293 


-1 


.91359 


0, 


.00860 


-0.02245 


4945 


.799895 


-1 


.85815 


0, 


.00807 


0.01502 


4947 


.745317 


-1 


.93960 





.00741 


-0.03134 


4948, 


.673112 


-1 


.82992 


0, 


.00743 


-0.06231 


4972, 


,707323 


-1 


.93123 


0, 


.00854 


-0.03631 


4975. 


,733486 


-1 


.93144 


0, 


,01100 


-0.01416 


4982, 


.647535 


-1 


.83433 





,01036 


-0.02677 


4984. 


.642389 


-2 


.04210 





,00892 


-0.04270 


4985. 


.694776 


-1 


.81561 


0, 


,00802 


-0.02406 
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The analysis was perfor med in a very s imilar 



fashion to that described bv lWest et al.l (J2009|) us- 
ing a spectral synthesis package and LTE model 
atmospheres. The H a and Hp lines were used to 
determine the effective temperature (T e g), while 
the Na I D and Mg I b lines were used as surface 
gravity (log 3) diagnostics. In addition the Ca H 
& K lines provided a further check on the derived 
T e g and \ogg. The elemental abundances of sev- 
eral elements were determined from measurements 
of several clean and unblended lines. The param- 
eters and abundances obtained from the analysis 
are listed in Table [2] 

In our spectra the Li 1 6708A line is not detected 
(EW < 2mA) , allowing us to derive an upper-limit 
on the Lithium abundance of log n(Li/H) + 12 
< 0.8. The lack of lithium would imply a n age 
in excess of 5Gyr (jSestito fc Randichl 120051 ). The 
stellar rotation velocity (v sin i) was determined 
by fitting the profiles of several Fe 1 lines using an 
average value of u mac = 2.0 kms -1 for the macro- 
turbulence (w mac )- 

In addition to the spectral analysis, we have 
also used available broad-band photometry to es- 
timate the total obser ved bolometric flux. The 
Infrared Flux Method (|Blackwell k Shallisl li~977r i 
was then used with 2MASS magnitudes to deter- 
mine T c ff and stellar angular diameter (6). This 
gives T e ff= 5550 ± 130 K, which is in close agree- 
ment with that obtained from the spectroscopic 
analysis (T off = 5700 ± 150 K). 

Co mparison with the s tellar evolution mod- 
els of iGirardi et al. (|2000l ) for solar metallicity 
(Z = 0.02) gives maximum-likelihood values A/* = 
I.OOIq^M q as shown in Figur e H Alternative 
models from iBaraffe et aT ( 19981) give essentially 
the same results as the stellar evolution models 
have close agreement in this mass range. The un- 
certainties on the stellar density le ad to a large 



uncer tainty on the age from the IGirardi et al 



(2000) isochrones producing an estimated age of 



2.3 



+5.8 



Gyr. 



3.2. Planet parameters 

The CORALIE spectroscopic RV data were 
combined with the WASP-South and EULER- 
CAM photometric data in a simultaneous fit to 
determine the planetary parameters. The method 
of Markov Chain Monte Carlo (MCMC) as de- 



tailed in previous investigations (|Pollacco et al 
20081 : ICollier Cameron et alJ hoOlh was used. We 
use the Claretl ( 2000 ) limb darkening coefficients 
for the appropriate stellar temperature and pho- 
tometric passband and a adaptive stepsize mech- 
anism is used during the 5000 step burn-in phase 
until the chain converges. At the end of this phase, 
the adaptive stepsize mechanism is switched off for 
the final 20000 steps in the chain. The autocor- 
relation length of the chain was 9 ± 1 for all the 
parameters, indicating that no unwanted correla- 
tions are present and the chain is "well-mixed" . 

Initial fits showed that the eccentricity was 
poorly constrained but consistent with zero and 
so was fixed at this value in subsequent fits. The 
prior on the stellar mass was set to 1.0 M©, as 
indicated by the evolutionary tracks discussed in 
the previous section, but no constraint or prior on 
the stellar radius or density was used in the fit. 
The transit parameters such as the period, depth, 
duration were initially set at the values from the 
transit search of the WASP-South data and sub- 
sequently refined in the MCMC code using all the 
available data. 

The best fitting system parameters are listed in 
Table [3] and show that WASP-16b is a reasonably 
close Jupiter analog albeit somewhat less massive 
and in a P ~ 3 day orbit. The host star has a 
fitted mass and radius which are slightly smaller 
than the Sun, leading to a slightly higher density 
than the solar case but all the parameters are iden- 
tical to the Sun within the error bars. The lack 
of lithium detection, low v sin i and similar large 
inferred age also point towards WASP- 16 being a 
solar analog hosting a hot Jupiter planet. 

4. Times of transit 

Although WASP-16 was observed with WASP- 
South for one full and two partial seasons, there 
are very few complete transits within the time- 
series suitable for determining times of transits. 
This illustrates the need for long t imeseries on 



poten tial transit fields as shown by iSmith et al 



(|2006h . In total we find four complete and well 
measured transits from the SuperWASP data and 
these are shown in Table 2] as 'Fitted T a ' along 
with one time of transit determined from the EU- 
LERCAM data. The predicted times of transit 
from the MCMC ephemeris (given in Table [3]) are 
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Table 2 

Stellar parameters of the WASP- 16 host star. 



Parameter Value 
R.A. = 14 h 18 m 43!92, Dec = -20°i6'31"8 (J2000.0) 





5700 ± 150 K 




4.5 ± 0.2 




1.1 ± 0.2 kms" 


v sinz 


3.0 ± 1.0 kms" 


Spectral Type 


G3V a 


[Fe/H] 


0.01 ± 0.10 


[Na/H] 


0.15 ± 0.08 


[Mg/H] 


0.14 ± 0.10 


[Si/H] 


0.10 ± 0.07 


[Ca/H] 


0.11 ± 0.12 


[Sc/H] 


0.14 ± 0.07 


[Ti/H] 


0.05 ± 0.14 


[V/H] 


0.09 ± 0.15 


[Cr/H] 


0.02 ± 0.11 


[Co/H] 


0.17 ± 0.08 


[Ni/H] 


0.07 ± 0.12 


logA(Li) 


< 0.8 


T eff (IRFM) 


5550 ± 130 K 


6»(IRFM) 


0.052 ± 0.003 m 



a Estimated from J-H color 
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Table 3 

System parameters for WASP- 16b. 



Parameter 



Value 



Error 



P (days) 
T (HJD) 

Tdur (days) 

i?p/R* 

b = acosi/R* 



3.1186009 
2454584.42878 
0.0800 
0.01199 
0.798 



+0.0000146 
-0.0000131 
+0.00035 
-0.00023 
+0.0018 
-0.0012 
+0.00052 
-0.00039 
+0.026 
-0.019 



K x (km s- 1 ) 
7 (km s _1 ) 
a (AU) 
i (degs) 



(adopted) 

1167 +0.0024 
u - 110 ' -0.0019 

-1 93619 +0.00021 

i.aouia -0.00023 

0421 +0.0010 

-0.0018 

85.22 ±Ul 



M* (M ) 
R* (R©) 
log g* (cgs) 
p* (p©) 



1.022 
0.946 
4.495 
1.21 



+0.074 
-0.129 
+0.057 
-0.052 
+0.030 
-0.054 
+0.13 
-0.18 



M p (M Jup ) 0.855 

R P (Rjup) 1-008 

P P (pjup) 0.83 

log 5p (cgs) 3.284 

T eq (A = 0,F = 1) (K) 1280 

Safronov number (6) 0.070 



+0.043 
-0.076 
+0.083 
-0.060 
+0.13 
-0.17 
+0.041 
-0.064 

+35 

-21 

±0.010 
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also shown in Table [4j There is currently an in- 
sufficient number of measured transits with ade- 
quate precision to suggest anything other than a 
constant period. 

5. Conclusions 

We report the discovery of a new transiting 
planet with the WASP-South station of the Super- 
WASP survey. The planet, designated WASP-16b, 
orbits a star which is a close solar analog, having 
temperature, mass, radius, metallicity and grav- 
ity the same as the Sun, within the error bounds. 
The age of the host star, with an admittedly large 
error bar, is also close to the solar age. 

The orbiting planet is a reasonable Jupiter ana- 
log although somewhat less massive than Jupiter 
(M p ~ 0.85 Mj up ), but with a near identical ra- 
dius (R p ~ 1.01 Rj U p) leading to a density some 
80% of Jupiter. This planet falls in the lower 
left corner of the group of "normal" Jupiter-sized 
planets in the Mass/Radius diagram, with the 
majority of objects in this region being some- 
what larger than Jupiter in either mass or ra- 
dius. Additionally if we compute the Safronov 
number 6 = \{V esc /V orb ) 2 = 0.070 ± 0.010 
for this planet, this along with it's equilib- 
rium temperature T eq — 1280 K places it in 
the center of the Class I planets as defined by 
Hansen fe Barmanl ( 2007 ). The "normality" of 
this planet makes it similar to WASP-2b, TrES-lb 
and other "normal" extrasolar planets ans stands 
in contrast to the inflated radii and low dens i- 
ties of planets like TrES-4b (ISozzetti et al.ll2009f) . 
HP 209458b dBrown etUI l200lh and WASP-lb 
(ICollier Cameron et al.ll2007l ). 
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ties of Keele, Leicester, St. Andrews, the Queen's 
University Belfast, the Open University and the 
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0.9 0.925 0.95 0.975 1 1.025 1.05 1.075 1.1 
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Fig. 1. — Zoom of the transit region of the Su- 
perWASP light curve of WASP-16b with the best 
fitting MCMC model overplotted 
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Fig. 3. — Radial velocity curve (upper panel) of 
WASP-16 from the Swiss 1.2m+CORALIE along 
with the best-fitting model which includes the pre- 
dicted Rossiter-McLaughlin effect. The resulting 
bisector spans are shown in the lower panel. The 
uncertainties on the bisector spans are double the 
radial velocity uncertainties. 




0.99 1 1.01 1.02 1.03 1.04 1.05 
Transit phase 



Fig. 2.— I c band light curve from EULERCAM 
and residuals from the transit fit of WASP-16b 





Fig. 4.— The posi tion of the WASP-16 host star 
in the isochro n es of l Girardi et alj (|2000h (left) and 
Baraffe et ail (|l998| ) (right). In both plots the 
0.9, 1.0 and 1.1 M©mass tracks are shown along 
with 100 Myr (solid), 1 Gyr (dashed), 5 Gyr (dot- 
dashed) and 10 Gyr (dotted) isochrones. 
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